Finite Size Scaling for the Glass Transition: the Role of a Static Length Scale 
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Over the last decade computer simulations have had an increasing role in shedding light on difficult 
statistical physical phenomena and in particular on the ubiquitous problem of the glass transition. 
Here in a wide variety of materials the viscosity of a super-cooled liquid increases by many orders 
of magnitude upon decreasing the temperature over a modest range. A natural concern in these 
computer simulation is the very small size of the simulated systems compared to experimental ones, 
raising the issue of how to assess the thermodynamic limit. Here we offer a theory for the glass 
transition based on finite size scaling, a method that was found very useful in the context of critical 
phenomena and other interesting problems. As is known, the construction of such a theory rests 
crucially on the existence of a growing static length scale upon decreasing the temperature. We 
demonstrate that the static length scale that was discovered in Ref. [3] fits the bill extremely well, 
allowing us to provide a finite size scaling theory for the a relaxation time of the glass transition, 
including predictions for the thermodynamic limit based on simulations in small systems. 
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I. INTRODUCTION 

Glasses are ubiquitous in nature; the problem of the 
glass transition where the viscosity of a super-cooled liq- 
uid increases by many orders of magnitude upon a tem- 
perature decrease of a couple of hundred degrees remains 
one of the major challenges in modern condensed matter 
physics. With the advent of powerful computers, numeri- 
cal simulations played a vital role in revealing some of the 
crucial characteristic associated with this phenomenon 
0, dj. In spite of this progress, due to still present lim- 
itations in computer power, simulations are still done in 
small sample sizes from thousands to a few tens of thou- 
sands of particles. Compared to experimental samples 
this is minute. Thus finite size effects in the results are 
always a major concern, especially when one deals with 
a system where a characteristic length scale grows with 
decreasing temperature p|. 

Finite size effects are not always unwanted; in fact one 
can extract useful insights by studying these effects care- 
fully [4|. As an example the temperature dependence of 
the static length scale can be obtained from the system 
size dependence of some relevant observable by analyz- 
ing the data using finite size scaling. In the glass com- 
munity the existence and usefulness of a static length 
scale is still debated; even if one can extract some static 
length scale by doing some sophisticated analysis 0-0]; 
the slow growth of this length scale makes it difficult to 
extract reliable insights . Nevertheless finite size effects in 
the dynamics of the supercooled liquids certainly exist, 
notwithstanding the fact that the characteristic length 
scale may not change very much. A careful analysis of 
these effects can shed a substantial light on the physics 
of the glass transition [6l4l9||. 

In Q, we already showed that our proposed length 
scale is the same object which governs the finite size ef- 
fects exhibited in configurational entropy S c in the Kob- 
Andersen Model in 3 dimensions [7[ . Further more t q and 
S c are related via the empirical Adam-Gibbs relation [2(| 



or may be via more sophisticated RFOT Theory [2lM23j . 
So our length scale directly relates dynamics with ther- 
modynamics in glass forming liquids via the Adam-Gibbs 
relation. Unfortunately we can not determine the expo- 
nents predicted in the RFOT Theory since the range of 
increase of our length scale is not sufficient to determine 
reliably the exponents. 

In this study our main aim is to establish that our 
proposed length scale is a crucial object for the theory 
of glass transition by demonstrating that it faithfully ex- 
plains the finite size effects seen in relaxation times. In a 
similar attempt in [l6[ to explain the finite size effect seen 
in the relaxation time, it was argued that a large system 
of particles can be thought of as a collection of smaller 
subsystems called "elementary units" which are coupled 
to each other by coupling constants. The finite size ef- 
fects seen in the relaxation times can be thought of as a 
organized rearrangements of these elementary units. It 
is quite possible that our proposed length scale actually 
measures the size of these elementary units, but further 
study is needed to pinpoint this connection. It is also 
worth mentioning that in [l6j a linear relationship has 
been proposed between the relaxation time and the sys- 
tem size for small system sizes. In contrast we propose a 
scaling function associated with the underlying growing 
static length scale. 

A clear advantage of our proposed length scale Q is 
that it is relatively easy to extract both in simulations 
and in experiments |24|_a,s compared to the other length 
scales like point-to-set Q and patch length scale [25|, [2|| . 
While we can compute our length scale for larger range of 
temperature this is prohibitive for the other length scales, 
thus at this time there is not enough data to provide a 
conclusive comparison between these length scales. 

Below we present a systematic study of the system size 
dependence of the a-relaxation time of different model 
glass forming systems in both two and three dimen- 
sions. We choose two completely different interaction 
pair potential- one with pure repulsion and one with 
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short range repulsion and long range attraction to stress 
the generality of the finite size effects in the dynamics of 
glass forming systems 



II. MODELS AND METHODS 

Here we study two distinct models, one with pure 
power law repulsive interaction d, [28| and the second 
being the well-studied Kob- Andersen Model [2!| in both 
2 and 3 dimensions. We will denote the Kob-Andersen 
Model in 3-dimensions as 3d KA and the slightly modi- 
fied version of it for 2-dimensions as 2d mKA [30]. The 
pure repulsive models will be referred to as 3d RIO and 2d 
RIO for 3-dimensions and 2-dimensions respectively. The 
usual Kob-Andersen model with a 80 : 20 binary mixture 
in 2-dimension shows prominent clustering effects so the 
model was slightly modified and a 65 : 35 binary mix- 
ture was used. This modified model shows no clustering 
behavior in all the simulated temperature range [30l l3l| . 
Throughout this paper we use temperature units such 
that the Boltzmann constant equals unity. 

The potential of the Kob-Andersen model is given by 



(1) 



where a,/3 G {A,B} and caa = 1.0, £ab = 1-5, (bb — 
0.5, <r AA = 1-0, a A B = 0.80, obb = 0.88. The Interac- 
tion Potential was cut off at 2.50ct Q( 3. We have performed 
the simulations at six different temperatures in the range 
T G [0.60, 0.43] in 3-dimensions and T G [0.90, 0.45] in 
2-dimensions at a number density p — N/V = 1.20. 

The interaction potential for the pure repulsive model 
is given by 
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where ry is the distance between particle i and j, e is the 
energy scale, and x c is the dimensionless length for which 
the potential vanishes continuously with q derivatives. 
The interaction length scale Ay between any two particles 
i and j is Ay = 1.0 A, Ay = 1.18 A and Ay = 1.4A for two 
'small' particles, one 'large' and one 'small' particle and 
two 'large' particle respectively. The coefficients cut are 
given by 
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We chose the parameters x c = 1.3854, k — 10 and q = 2. 
The unit of length A is set to be the interaction length 
scale of two small particles, and e is the unit of energy. 
We have performed the simulations at six different tem- 
peratures in range T G [0.80, 0.52] at number density 
p = 0.81 in 3-dimensions and T G [0.65, 0.48] at p = 0.85 
for 2-dimensions. 



We have simulated these models in the NVT ensem- 
ble using a Berendsen thermostat [32| to maintain the 
temperature. The system size was in the range N G 
[128, 20164] and we used as many as 20 to 50 different 
samples to average the data depending on the system 
size. Generally for small system sizes and at low tem- 
perature the fluctuations increase, demanding more ex- 
tensive averaging of the data to get reliable estimates of 
the wanted observables. For each system size and each 
temperature we measured the a-relaxation time and the 
characteristic length scale. The former is calculated from 
the decay of the overlap correlation function defined as 
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where the weight function w(x) = 1 if x < 0.30 and zero 
otherwise. The a-relaxation time r Q is the time where 
this correlation function goes to 1/e of its initial value i.e. 
Qija) = l/ e - The characteristic length scale was mea- 
sured using the method proposed in Ref. [3j . The reader 
is encouraged to examine this reference for the details 
of the method. In brief, it is based on the insight that 
the minimum eigenvalue of the Hessian matrix belongs 
to either an elastic mode or to a plastic mode, depend- 
ing on the system size. The characteristic length is the 
crossover system size between these two possibilities. The 
Hessian matrix is obtained for every supercooled system 
at temperature T by employing the conjugate gradient 
minimization methods to get an inherent structure (at 
temperature T — 0), and then using the Lanczos method 
(33j to get the minimum eigenvalue. 



III. RESULTS 

In the left panel of Fig[TJ we present the system size 
dependence of the a-relaxation time r a for the Kob- 
Andersen model in 3-dimensions, for different temper- 
atures. The data indicates a substantial system size de- 
pendence of the relaxation time which appears to increase 
slowly with decreasing the temperature. This is an indi- 
cation that the length scale which governs this behavior 
increases also with decreasing the temperature. Phys- 
ically we expect the relaxation time to decrease with 
increasing system size if one asserts the existence of a 
cooperative length scale. This length scale becomes rel- 
evant to the system's dynamics when the temperature 
goes bellow the so-called 'onset temperature' j35( which 
separates the characteristic Arrhenius regime from the 
'fragile' regime of the temperature dependence of the re- 
laxation time, cf. 34|. The argument is as follows : For 



a system of size N at temperature T, the typical free 
energy barrier it needs to cross for relaxation is 

AF(N,T) = AE(T)-Tlog[g(N,T)}, (5) 

where AE(T) is the potential energy barrier. We ex- 
pect the potential barrier to be determined by the en- 
ergy landscape which in general does not depend on the 
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FIG. 1: (Color online) 3d KA Model : Left Panel: The system size dependence of the a-relaxation time for the Kob- Andersen 
Model in 3-dimensions for different temperatures. The lines are fits to the data using the fitting function f(x) — A + B/x to 
extract the large N value T a (N — > oo,T) of the relaxation time. Middle Panel: The y-axis is rescaled by this value such that 
all the data for large N converge to unity. Note that this procedure does not collapse the data indicating the existence of a 
length scale. Right Panel: The complete collapse of the data using the length scale obtained in 0] (see text for details). 



system size. On the other hand \og[g(N,T)} is the en- 
tropic contribution to the free energy barrier. Now at 
temperatures lower than the onset temperature [35| the 
cooperative relaxation mechanism plays the main role for 
relaxation and one expects the degeneracy factor for this 
case to be system size dependent when the system size is 
small. 

Consider a system of N particles in a cubic box of size 
L at some temperature T. In the non-Arrhenius regime 
with a typical length £(T), we expect for a system with 
L << £(T), that any relaxation event will occur on the 
scale L in order to be successful. So the degeneracy factor 
will grow with system size. For a system with L » £(T) 
we expect g(N, T) ~ const without any change when the 
system size grows. Thus the free energy barrier can be 
represented as a scaling function of L /£ (T) according to 



AF(N,T) = AE(T) - Tlog 
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where A is a proportionality constant independent of the 
temperature, and obviously 



f(x) — > Const(T) when x — > oo 



(7) 



Next we calculate the system size dependence of the re- 
laxation time for different temperatures, using the es- 
timate r a « To expJAi 7 '(iV, T)/T] where To is a typical 
attempt time, cf. [34| . The result is again a scaling func- 
tion of L/£(T), reading 
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From the property ([7]) of the scaling function we con- 
firm that t(N,T) -> ConstlT) when N -> oo. 

Examine the consequences of these arguments for the 
data of the Kob-Andersen model, cf Fig. [TJ In the 
left panel we present the a relaxation time for differ- 
ent temperatures and systems sizes. In the middle panel 
of Fig[T] we rescaled the y-axis by the large N value of 
of the relaxation time r a (oo,T) for the different tem- 
peratures. We estimated this value by fitting the data 
with a functional form f(x) = a + b/x ; One sees that 
this fitting is adequate for obtaining the large TV value 
of the relaxation time. In the right panel of FigQ] we 
present the full data collapse which is obtained by reseat- 
ing the system size using the length scale obtained in Q . 
We stress that no adjustments were necessary and none 
were made to these length scales, they were employed 
as obtained. Now one may ask whether a similar col- 
lapse can be done using the other static length scales 
mentioned above. In a recent study [36[ it was shown 
that the point-to-set length scale in the Kob-Andersen 
model in three dimension changes by a factor of 2.2 in 
the temperature range T 6 [1.00,0.550]. Assuming that 
we can extrapolate these data to lower temperature we 
find point-to-set length scale changes by similar factor 
in the temperature range of our study. Thus it appears 
possible to collapse the relaxation time data using the 
point-to-set length scale in this temperature range. This 
suggests that point-to-set length scale is proportional to 
our length scale. It remains a future work to compute 
point-to-set length scale to lower temperature to confirm 
this connection. 

In Fig 12] we present the results of a similar analysis 
for the model system with pure repulsive interaction (3d 
RIO) in 3-dimensions. The data collapse observed in the 
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FIG. 2: (Color online) 3d RIO Model: Left Panel: The system size dependence of the a-relaxation time for a system interacting 
via a 1/r 10 pure power law repulsive potential in 3-dimensions for different temperatures. The lines are fits to the data using 
the fitting function f(x) — a + b/x to extract the large N value T a (N — > oo,T) of the relaxation time. Middle Panel: same 
rescaling as in the previous figure. Right Panel: The complete collapse of the data using the length scale obtained in Q (see 
text for details). 
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FIG. 3: (Color online) 2d R10 Model: Left Panel: The system size dependence of the a-relaxation times for a system interacting 
via a 1/r 10 pure power law repulsive potential 2-dimensions for different temperatures. The middle and right panels repeat the 
procedure explained in Fig. 1 











■ 


A 


o 








A 





o 






A 


A 


V 






10 



10 



N 



o 

A 




10 10 




v- 



• 


T = 


0.450 


□ 


T = 


0.500 





T = 


0.550 


A 


T = 


0.600 


V 


T = 


0.700 


< 


T = 


0.900 



10' 10° 

N/£ d (T) 



^0 A 



FIG. 4: (Color online) 2d mKA Model: Left Panel: The system size dependence of the a-relaxation times for modified Kob- 
Andersen Model (mKA Model, see text for details) in 2-dimensions for different temperatures. The fits and the middle and 
right panels follow the procedure described in Fig. 1 
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FIG. 5: (Color online) Temperature dependence of the static 
length scale £ for different models studied. 
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FIG. 6: (Color online) Relaxation time in the thermodynamic 
limit plotted according to Eq. QT] ( see text for details) for 
all the four studied models. The apparent universality is en- 
couraging. The parameter Ao used here (see text for details) 
is model dependent and numbers are quoted in the text. 



right panel is again quite good. We present similar anal- 
ysis in 2-dimensions for the model with pure repulsive 
interaction in Figj3] and for the modified Kob- Andersen 
model in FigfJ] In FigfSJ we plotted the temperature 
dependence of the correlation length for all the models 
studied to compare the model to model variation of the 
growth of this length scale with decreasing temperature. 



IV. UNIVERSAL RELATION BETWEEN r a 
AND £ 

Having found that the typical scale fits the bill for finite 
size scaling, we should ask next whether there is a direct, 
and possibly universal relation, between the length scale 
and the time scale. We start with the following ansatz for 
the relation between structural relaxation time r„ in the 



thermodynamic limit and the static correlation length £ 

'Ao^(T)^ 



r a (T) oc exp 



T 



(9) 



where Ao is a non-universal coefficient that depends on 
the details of the glass former. To set a scale to this 
relation we choose a reference temperature T — Tq where 
the typical length is £o = 1.0. Then the relaxation time 
at that temperature is 



T a (T ) oc exp 
So we have the following relation 
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We reiterate that in the above equation the pre-factor 
s = An i s no t known apriori for different models. In 
Fig|5] we plotted the relaxation time data according to 
EqUU with A = 1.0, 0.70, 1.0,and 0.4274 for KA model 
in 3D, R10 model in 3D, R10 model in 2D and modified 
KA model in 2D respectively. We choose -0 = 1-0 for 
all the models. The collapse of the data indicates an 
encouraging possible universal relation between the static 
length scale and the relaxation time. 

To conclude, we have studied the finite size effects in 
the dynamics of supercooled liquids for different models 
systems in both 2 and 3 dimensions. In this article we 
have focused on the system size dependence of the a- 
relaxation time and found that the finite size effects can 
be very well explained by the static length scale that we 
proposed in Q for all these different models. In a recent 
study [13], it was argued that a growing static length 
scale is necessary to explain the finite size effects seen in 
the simulational studies of model glass forming liquids. 
We want to stress that our proposed length scale @ is 
indeed that length scale which faithfully explains all the 
finite size effects seen in the simulation results. 

Further more we realized that there is a possibility of 
universal relation between this length scale and the a- 
relaxation time when relaxation time for different models 
are plotted as a function of the length scales, all the data 
fall nicely on a master curve as also recently reported in 
[3^ | . This connection again seems to suggest that point- 
to-set length scale might be proportional to our proposed 
length scale. We reiterate that further study is needed to 
confirm this relation and care should be taken as these 
length scales and the relaxation times are not varying by 
much in the temperature range for which we could do the 
simulations. We hope that our work will inspire other re- 
searchers to do further studies to delineate, especially in 
the low temperature range, the validity of this apparent 
universal behaviour between the relaxation time and the 
static length scale. 
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